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Abstract 

By solving the two-body Schordinger equation in a very high precise nonperturbative numerical 
(NPnum) way, we reexamine the contributions of fine, hyperfine structure splittings of muonic 
hydrogen based on the Breit potential. The comparison of our results with those by the first order 
perturbative theory (^®*PT) in the literature shows, when the structure of proton is considered, 
the differences between the results by the and NPnum methods are small for the hne and 

hyperfine splitting of 2P state, while are about —0.009 meV and 0.08 meV for the F = 1 and 
total hyperfine splitting of 25 state of muonic hydrogen, respectively. These differences are larger 
than the current experimental precision and would be significant to be considered in the theoretical 
calculation. 
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I. INTRODUCTION 


Q 


In 2010, a precision measurement [l| of the Lamb shift in muonic hydrogen by using pulsed 
laser spectroscopy was performed and gave = 206.2949 meV. Combing this 

precise value with the theoretical calculation 


r 


■' 1/2 


= 209.9779 - 5.2262r2 + 0.0347r3 

^-^ 3/2 ^^ 1/2 ^ ^ 


(1) 


the values of the proton radius is extracted as Vp = 0.84184 fm [l|. In 2013, the further 
precise measurements of 2S — 2P transition frequencies of muonic hydrogen gave the 
magnetic radius of proton tm = 0.87 fm and the charge radius te = 0.84087 fm which are 
not signihcantly different from the value given by Ref. [l|. 

On the other hand, based on the hydrogen data or the ep scattering data, CODATA- 
2010 gave Vp ^ 0.878 fm j^, which is much larger than the results by the muonic hydrogen’s 
Lamb-shift. And if this value of proton radius is used, the theoretical prediction for the 
Lambs shift of muonic hydrogen gives jd] 


-l^npF=2 

3/2 


-E'^F=i|rp=o. 878 /m — 205.9726 meV, 


( 2 ) 


which deviates from the experimental Lamb shift of muonic hypdrogen about 0.32 meV. 

Many theoretical calculations j^, data analysis jo] and possible new mechanisms such as 
the three body physics jjl, the new exotic particles interactions js], the higher-order contri¬ 
bution of the hnite size (9| etc., have been discussed to try to understand such discrepancy. 
And also new experiment of ep scattering is proposed in JLab Ij. Combining all these 


current analysis, briefly, the radius of proton is still not well understood. 

For the muonic hydrogen, the energy transition of 2Py2‘^ and 2S^i^^ usually are expressed 
as 


T?Th 

■^OPF=2 

^-^3/2 


= ae 


2S-2P 

Lamb 


AE% + 


^HFS 




( 3 ) 


In the literature, the contributions of the four terms are usually calculated by the per¬ 
turbative theory. Using the quasipotential method in quantum electrodynamics ll|, the 
contributions to the four terms can be expressed as QQQ 
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FIG. 1; The 2S and 2P energy levels of muonic hydrogen. 
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where AEovp (205.0074 meV by the first order perturbative theory (^®*PT) and 0.1509 meV 
by the second order perturbative theory), AEks (1.5081 meV) and AEse (—0.6677 meV) are 
the energy shifts due to the one-loop vacuum polarization, two-loop vacuum polarization and 


the sum of self-energy and muonic-vacuum polarization, correspondingly, AEqed (0.0586 


meV) is the energy shift due to all further QED corrections, and the last two terms in 


are relevant radius-dependent contributions[l2|, AE^p^g (22.8054 meV), AE^pg 


2S,B 


2^3/21^ 


(3.392112 meV) and AE'^g^ (8.329150 meV) are the Fermi energies, AE"^^^^ (0.0266 
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meV), (—0.000886 meV) and (0.017637meV) are the contributions 

from the anomalous magnetic moment of muon, AE’^^p^^'' , and AE^pg^^'°^^^'^ 

are the other contributions Q,Q. 

Some of the above perturbative results have been checked by the nonperturbative nu¬ 
merical (NPnum) calculations, for example within the framework of the multiconfiguration 
Dirac-Fock (MCDF) method in and shotting-like method using quad-precision Fortran 


m 


16|. In this work, by using the Mathematica, we present another high precise NPnum 


calculations (much more precise than the quad-precision) on the energy shifts , 

and with considering the effects from the proton structure. And 

as a comparison, also the calculation of 5Eovp is presented. 


II. FORMULA AND NUMERICAL METHOD 

The leading order contribution to the fine and hyperfine structure of muonic hydrogen 
due to the proton structure is from the two photon exchange diagrams. Since there is IR 
divergence in these diagrams and such IR divergence is not dependent on the proton structure 
(only dependent on its charge), the one photon exchange diagram should be considered 
together in some way to cancel such IR divergence. This leads to the complexity in the 
discussion of the numerical calculation, so at present we take the effective potential from 
the one photon exchange diagram as an example to discuss the difference between the i®*PT 
calculation and the NPnum calculation. 



FIG. 2: One photon exchange Feynman diagram considering the form factors of proton. 

The Feynman diagram for one photon exchange with the proton structure is showed as 
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Fig. 121 where the vertex of 'jNN is taken as 


iFo 




( 5 ) 


with Fi = {Ge + tGm)/{^ + t),F2 = {Gm - Ge)/(1 - r), r = Q'^/ArUp, Ge,m the electro¬ 
magnetic form factors of proton, rup the mass of proton, and q = Ps—pi the four momentum 
of exchanged photon. The correction to the Coulomb potential from this diagram was dis¬ 
cussed in the recent work 17|, and in this work, we discuss its corrections to the hne and 


hypehne Breit potential by a precise numerical method. By using the quasipotential method 
the Breit potential can be expressed as [^, f^ . 
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where the One-loop Uehling potential Vovp is also presented for comparison and 


( 6 ) 


Gfs = -(1 + - (1 + r^)T 


+ 




(1-A;2)2' 

P ^—mkr 


,^2g-mr 


(1 - k^y 




2 9 

hi = + mr) -e”™"- 


9 9 T T 

h 2 = + mr) - 


2 2 
m r 


hs = + mr) - 


hi — (1 T 2Kp)(l -\~ h^) -\~ 


Hr. 


(1-F)2 


e-’""(l + mr) + 


Kp 


Hr 


1-k'^ 2 


(1-F)2 


e-”^'="(l + mA:r), 


5 




























3(cr^ • r)(crp • f) 


2L • (Tp 

9 qF=0 
^'^1/2 

-3 

-3 

0 

ooF=l 

^'^1/2 

1 

1 

0 

9 pF=0 

^W/2 

-3 

1 

-4 

9 pF=l 

1 

1 

3 

4 

3 

2pF=2 

3 

5 

1 

2 

9pF=l 

^^3/2 

-1 

5 

3 

10 

3 


TABLE I: The contributions of the spin related operators to 2S and 2P states. 


with a the fine-structure constant, jip and /i^ the anomalous magnetic moments of proton 
and muon, m the parameter in the electromagnetic form factors of proton which can be 
related with proton size as = 12/m^, k = 2mp/m, me, m^ the mass of electron and 
muon. To include the effects from the proton structure in the above Breit potential, the 
electromagnetic form factors of proton Ge,m{Q'^) have been taken approximately as the 
usual dipole form as [18 1 


— GM(Q^)/hp — 1/(1 + 


2\2 


(7) 


with = —q^. The contributions of spin related operators in the Breit potential to the 2S 
and 2P states are listed in Tab. HI 

When the effects from the proton structure are neglected in the above effective potentials 
(taking m —)■ oo), the energy shifts by i**PT reproduce the same 


and with those used in the literatures 


3,Q. 


The corrections from the proton 


structure in the above Breit potential is corresponding to replace the zero momentum transfer 
approximation by including the full dependence of the electromagnetic form factors of 
proton in the one photon exchange Feynman diagram. This is different with the corrections 
from the two photon exchange diagrams discussed in [^, 141. Since in this work our focus 
is on the difference between the results by i®*PT and NPnum calculation, we take the Breit 
potential as an example to show the difference. 

In our numerical calculation, we use the shotting method to hnd out the energy spec¬ 
trum by solving the reduced Schrodinger equations for u{r) directly, with u{r) = R{r)r, 
'0(r) = R{r)Y{6, 0) and '0(r) the wave function. In the detail, for muonic hydrogen we take 
approximately u{r = 0) ~ u{r = and u{r = oo) ~ u{r = 10®/m) to simulate 
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the behaviors of wave function at the boundary. We keep 200 digits of the numbers in the 
calculation and take the PrecisionGoal of NDSolve in Mathematica as 15, respectively. By 
these approximations, as a check we reproduce the energy spectrum of muonic hydrogen 
under the Coulomb potential with the precision better than meV. Since the numerical 
calculation is based on the shotting method, the precision is not sensitive on the form of the 
potentials or the solutions, this is different with the basis expansion method or variational 
methods usually used, and ensures our numerical calculation reliable for the other potentials. 


III. NUMERICAL RESULT 


In Tab. U and m we present the results by the and NPnum calculation including 
the effects from the proton structure in the Breit potential. The corresponding results 
without considering the proton structure by i®*PT and from the one-loop Uehling potential 
are also presented for comparison. 


rp(fm) 

AEjTs 

A p2P3/2 

AE 

^■^ovp 


NPnum 

Istprp 

NPnum 

Isiprp 

NPnum 

0 

8.34676 

- 

3.39121 

- 

205.00659 

205.15747 

0.83112 

8.34673 

8.34703 

3.39121 

3.39123 

205.00659 

205.15747 

0.84184 

8.34672 

8.34702 

3.39121 

3.39123 

205.00659 

205.15747 

0.87800 

8.34672 

8.34702 

3.39121 

3.39123 

205.00659 

205.15747 


TABLE 11: Energy shifts of different potentials using perturbative and precise numerical calcula¬ 
tions in meV where ^^*PT and NPnum denote the first order pertubative and precise nonperturba- 
tive numerical calculation, respectively. The typical proton size = 0,0.83112.0.84184,0.878 fm 
are taken as examples for comparison. The results by ^^^PT are same with those in Q,Q. 


From the last two columns of Tab. m we see the results by i^*PT and NPnum calcu¬ 
lations give about 0.15088 meV difference for Vovp- Actually the second order perturbative 


calculation of W™ gives the contribution about 0.1509 meV 


12| . so the combination of the 


first and second order pertubative calculation of Vo^p is almost same with our NPnum cal¬ 
culation. When including the effects from the proton structure and taking rp = 0.83112 fm, 
0.84184 fm, 0.878 fm as examples, the differences between the i®*PT and precise NPnum 
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rp(fm) 

^Ejfpg{F — 1) 

^E‘jfpg{F — 0) 




NPnum 


NPnum 

Istprp 

NPnum 

lstpr£ 

NPnum 

0 

5.70798 

- 

17.12394 

- 

3 

- 

22.83192 

- 

0.83112 

5.67921 

5.66943 

17.03762 

17.12628 

3 

3.0208 

22.71683 

22.79571 

0.84184 

5.67884 

5.66918 

17.03651 

17.12406 

3 

3.0206 

22.71535 

22.79324 

0.87800 

5.67759 

5.66832 

17.03277 

17.11676 

3 

3.0197 

22.71036 

22.78508 


TABLE III: Energy shifts of different potentials using perturbative and precise numerical calcula¬ 
tions in meV. The notations are same with Tab. HI The results by i®*PT are same with those in 

Q. 


calculations are about 3 x 10“^ and 2 x 10“® meV for the fine and hyperfine splitting of 2P 
states, respectively, and also these two splittings are almost independent on the proton size 
in the region G [8.2, 9.0] fm. From the Tab. IIIIl we see the derivations of the and 

NPnum calculation are about —0.009 and 0.02 meV for the hyperhne splitting of 2S state 
AEjfpg{F = 1) and jAEjfpg, respectively, which should not be omitted comparing with the 
precision of current experiments. We want to emphasize that by the NPnum calculation, the 
ratios are not strictly equal to 3 as predicted by the ^^*PT, but are about 3.02 

as showed in Tab. im and the absolute difference between AEjfpg(E = 1) and jAE'jfpg are 
large. Such a discrepancy means the relation AEjfpg{E = 1) = jAEjfpg is not suitable to 
be used as usual. To show the results in a more direct way, we use the polynome of to £t 
the numerical results of AE]fpg(E = 1) in the region Vp G [8.2, 9.0] fm by taking one point 
every 0.001 fm and the results are expressed as 

AE'jfpsiE = 1, NPnum) = 5.68008 - 0.00261343rp - 0.0122699rJ meV 

AEfpsiE = 1,^"^ PT) = 5.70793 - 0.0346149rp + 0.0000613066rp meV (8) 

with the residual mean square of the fitting as small as about 10“^^ meV^. 

After including the difference of AE]fpg(E = 1) by our estimation, the theoretical energy 
shift Eq.([T]) is changed as 

E^PF=2 - = 209.9869 - b.22Q2rl + 0.0347r3 (9) 

and if r„ is taken as 0.878fm, then the modihed E'^A^^ — is estimated as 205.9816 

^^3/2 ^^1/2 

meV, which deviates from the experimental results about 0.31 meV. This means after using 


















the precise NPnum calculation of AEjfpg{F = 1), the discrepancy of the measurement of 
the proton size is reduced about 3%, while if we take the as input then the energy 

shift Eq.([I]) is changed as 


E^pF=2 - = 209.9579 - 5.2262r2 + O.OSdTr^ (10) 

and if Vp is taken as 0.878fm, then E'Fpp ^2 — E'FA^i is estimated as 205.9526 meV, which 

^^ 3/2 ^^ 1/2 

deviates from the experimental results about 0.34 meV. We see the discrepancy will be 
intensihed about 7%. The full results show we should be careful to deal with the splitting 
beyond the and should replace ^AEjfpg by AEjjpg{E = 1) in the calculation. 

According to Tab. [land Fig. [H we also have 


( 11 ) 


and 


AE — {E2pF=i E2gF=o) {E2pF=2 E2gF=i) — AE'jfpg AE^pg. (12) 

We note that AE depends on AEjfpg instead of AEjfpg(E = 1) appeared in E 2 pF^ 2 —E 2 gF=i, 
so for AE the NPnum calculation gives relative larger corrections. We present the results of 
AE from the Breit potential by the two methods in Tab. IIVI With the different proton radius 
as input, the differences of AE are about 10“^ ~ 10“^ meV. At present, the E 2 pF =2 —E 2 gF=i 
and E 2 pF=i —E 2 gF=o transition frequency in muonic hydrogen have been both measured with 
very high accuracy j^, the experimental value of AE is about 19.56 meV. Our numerical 
results are different from the experimental value since we have not considered the other 
corrections beside Breit potential. However, our calculation implies the precise NPnum 
calculation is needed. 

Since the obvious difference exists between the i^*PT and precise NPnum calculations of 
the hyperhne splitting of muonic hydrogen’s 2S state, we also present the similar comparison 
of the ep system in Tab. jV] where the results show the calculation for the hyperhne splitting 
of S state in ep system by the ^^^PT is also not good enough. Different with ^®*PT, the 
results by the precise NPnum method are more sensitive on the proton size. 

In a summary, by using shotting method in Mathematica we give a very high precise 
NPnum calculation for the energy shifts of the Breit potential including the effects from the 
proton structure. Our results show that when taking into account the proton structure, the 
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rp (fm) 

A£'(meV) 


NPnum 

0 

19.44070 

- 

0.83112 

19.32562 

19.40434 

0.84184 

19.32414 

19.40187 

0.87800 

19.31915 

19.39371 


TABLE IV: The difference of 2S — 2P transition energy AE of muonic hydrogen in meV with 
different radius of proton as input. The notations are same with Tab. HIl 


rp (fm) 




NPnum 

Istprp 

NPnum 

0 

1420478.8 

- 

177559.8 

- 

0.83112 

1420440.1 

1425769.6 

177555.0 

178221.2 

0.84184 

1420439.6 

1425700.8 

177555.0 

178212.6 

0.87800 

1420438.0 

1425481.1 

177554.7 

178185.1 


TABLE V: Frequencies in kHz of hyperfine splitting Ejfp^ and E'jfpg in hydrogen. 

precise NPnum calculations give very small corrections to the hyperfine splitting of 2 P 3/2 
and hne structure of 2P states, but give about —0.009 meV and 0.08 meV differences for 
the hyperhne splitting AEjfpg{F = 1) and AEjfpg of munoic hydrogen with those usually 
used in the literatures by The similar properties are also found in the hydrogen case. 
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